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Metric Problems for Quadrics in Multidimensional Space 
Alexei Yu. Uteshev0 and Marina V. Yashina^ 

St. Petersburg State University, St. Petersburg, Russia 

Abstract 

Given the equations of the first and the second order surfaces in M", our goal 
is to construct a univariate polynomial one of the zeros of which coincides with 
the square of the distance between these surfaces. To achieve this goal we employ 
Elimination Theory methods. The proposed approach is also extended for the case 
(NJ \ of parameter dependent surfaces. 

o_ 

^ ■ 1 Introduction 



We solve the problem of finding the distance d from the ellipsoid 

X^AiX + 2B'[X -1 = 



U: either to linear surface given by the system of equations 

C^X = 0,...,C^X = (2) 



or to quadric 

^ ■ X^AaX + 2B^X -1 = 0. (3) 

Here X = [xi, . . . , a;„]^ is the column of variables, {Bi, B2, Ci, . . . , Ck} C M" are the given 



CN I columns, Ai and A2 are the given symmetric matrices and Ai is sign-definite. 

^ ■ The distance is evaluated in Euclidean metrics || ■ II2, i.e. 

o 

(N 



d = min ^(X - - F) 

subject to {X e M"|X^AiX + 2BfX -1=0} and {Y G R"|Cf F = 0, . . . , C^Y = 0} or 
{Y E W\Y^A2Y + 2B^Y -1 = 0}. 

Being a problem of nonlinear optimization it can be solved via generation of a suitable 
iterative procedures [6], [8] or by application of some symbolic transformation of algebraic 
equations aiming at reducing the number of involved variables. Thus, for instance, for the 
distance problem between ([1]) and ([3]) the starting point is the following system resulted 
from the Lagrange multipliers method 

X-Y- Ai(AiX + Bi) = 0, -X + Y- X2{A2Y + B2) = O ... 
X^AiX + 2BfX = 1, r^AaF + 2B^Y = 1. 

Several publications [2], [6] were devoted to the development of mentioned approaches 
for the dimensions n = 2 and n = 3. They were focused on finding the coordinates X and 
Y of the nearest points in the considered surfaces. In comparison with those approaches, 
in the present paper we suggest an alternative one aimed first at evaluation of the distance 
itself. This is achieved via introduction of a new variable z by the equation 

z-{X- Yf{X -Y) = 0. 
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Being attached to the system (jl]) this equation provides the critical values of the distance 
function. For the obtained algebraic system one may apply an algebraic procedure of 
elimination of all the variables except for z. That means, it is possible to construct an 
algebraic univariate equation Ti^z) = one of the zeros of which (generically minimal 
positive) coincides with the square of the distance we are looking for. This construction 
can be performed either via the Grobner basis computation or with the aid of the classical 
Elimination Theory toolkit. We have chosen the second approach and succeeded in finding 
explicit expressions for the polynomial J-" for each of the stated problems in ra— dimensional 
space. Any of the real zeros of J^{z) = corresponds to a pair of points on the treated 
surfaces, and we also suggest an algorithm for evaluation of their coordinates. It turns 
out that these coordinates can be generically expressed as rational functions of the value 
of z. We also treat a surface intersection problem. Some of results from Sections [3] and H] 
were first formulated in [11] . In the present paper we give a proof for Theorem [6] (missed 
in pjj ) and correct one void in the proof of Theorem HI 



2 Algebraic preliminaries 

From the mentioned in the previous paragraph Elimination Theory toolkit, the most 
perfect gadget for our purpose turns out to be the discriminant. We will be in need of its 
univariate and bivariate form. 

Univariate discriminant. For the univariate polynomial g{x) = b^x^ + bix^~^ + 
. . . + e C[x], 6o 7^ 0, > 2 its discriminant is formally defined as 

N 

where {fii, . . . , /xat} is a set of zeros of g{x) counted in accordance with their multiplicities. 
We will also use an alternative definition of discriminant 

JV-l 

vAg) i-ir^^-^y^N^bt' n ^(^^O' (6) 

i=i 

where {Ai, . . . , Aat-i} is a set of zeros of g'{x) counted in accordance with their multiplic- 
ities. The constructive computation of discriminant - in the form of polynomial function 
of the coefficients of g{x) - can be performed with the aid of several determinantal repre- 
sentations. We will utilize the Bezout's approach [T] which is based on the coefficients of 
the remainders on dividing x^g{x) by g'{x): 

x^g{x) = beo + b^x + . . . + bi^N~2X^~'^ + qi{x)g'{x), qe{x) G C[x] 

for ^ G {0, . . . , — 2}. Compose the matrix from these coefficients 

® = [%]£-io- (7) 

Denote by ?Biv-i,j the cofactor to the corresponding entry of the last row of 
Theorem 1 One has 

P,.(^7) = ArX-'det®. 
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The polynomial g{x) possesses a multiple zero ijf det^ = 0. Under this condition, the 
multiple zero is unique iff 7^ 0; in this case it can be expressed rationally via the 

coefficients of g{x): 

A = (8) 



Example. Find the real values of the parameter a under which the polynomial 

g{x) =x^ + 6x'^ + 2x^ + ax'^ — x + 3 
possesses a multiple zero, and evaluate this zero. 

Solution. We compute first the remainders on division of g, xg, x'^g, x^g by g'{x): 

81 / 12a 4\ /3a 36\ . 124 o 
+ — — - t\ x+{— - — ]x^- — x-^, 



25 V 25 5 J V 5 257 25 

124 /248 a 81\ / 12 a 644\ . /3 a 2796\ , 
h \ x + { \ x^ + { \ x\ 

125 V 125 25/ V 25 125/ \ 5 125 J 

3 a 2796 / 6 , 5592 124\ /158 a 14751 \ n 

\ h a^ a ] x + { x^ 

25 625 V 25 625 125/ V 125 625 J 



84 a 63884 
+ ' "^5 62^' ^ ' 
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84 a 63884 f 168 , 128143 2796 \ 

Cfc H T^TTZ^ + ——— X 



125 3125 V125 3125 625 J 

f 6 2 3072 380204\ ^ /2nAa 1459461\ o 

+ a^ a H a; + \ x'^. 

V 25 625 3125 J \ 125 3125 J 

Then compose the matrix 53 from the coefficients of powers of x and compute its deter- 
minant 

, ^ (a + 7) (324 a^ + 5481 a^- 87771 a^- 409817 a + 5759315) 
® = 3125 • 

The discriminant 'Drc{g) coincides (up to a numerical factor) with the numerator of the 
last fraction and it vanishes iff 

a G {-24.63939477, -9.29644677, -7}. 

To evaluate the corresponding multiple zero of g{x), we utilize formula (jS]): 

27 , 18 2 32537 2724 

a H a H a H 

625 5 625 625 



54 4 1296 o 4508 2 17208 57532 
a a H a H a 



625 625 625 125 625 

where numerator and denominator are just the minors to the entries of the last row of 05. 
Substitution of the obtained values of a into this formula yields the corresponding values 
of muhiple zeros: -3.80947138, 0.74648466, -1. 

Corollary 1 Let (f){x) = p{x)/q{x) be rational function with relatively prime p{x) and 
q{x). Functions (f){x) and (j)'{x) posses a common zero iffVx{p) = 0. 
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Proof. One has = iff p{x) = 0. Let degp(x) = m and Ai, . . . , A„i stand for the 
zeros ofp{x). Thus 

,// N P'{x) p{x)q'{x) , P'i^i) r ■ n 1 

g(x) g^(a;) g(Aj) 
here g'(Aj) ^ under the assumption of the corollary. Therefore 

m mm 

j=i j=i j=i 

and in accordance with the definition (j6]) this product vanishes iff T>x{p) = 0. □ 

Corollary 2 For polynomial g{x) of degree N > 2 and a constant A G C one has: 

V,{A.g{x)) = A'^'-'VM, (9) 
V,{x-g{x)) = [g{0)]^VM- (10) 

Theorem 2 One can find polynomials providing the so-called linear representation of the 
discriminant, i.e., the pair {u{x), v{x)} C C[x] satisfying the identity 

v{x)g{x) +u{x)g'{x) = det^. (11) 

Here v{x) can be represented as the determinant of the matrix obtained on replacing the 
first column of 53 by [l,x, . . . , while 

where 5B denote the matrix obtained from 53 by replacing its first column by 

[0, bo,N-2, bo,N-2X + bi^N-2, bo,N-2x'^ + &l,Ar-2a; + 62,Af-2, • • • , 
bo,N-2X ^ + bi^N-2X + . . . + bN-'i,N-2\ ■ 

The polynomials u{x) and v{x) satisfy the restrictions 

deKu < N -1, dee:v < N - 2. 



Bivariate discriminant. For the given polynomial g{X) G C[X], X = {xi,X2), 
degg = N > 2 we define its discriminant as 

i=i 

Here Aj = [Xj^, Xj,^) £ stands for the stationary point of g{X), i.e. a zero of the system 
dg/dxi = 0, dg/dx2 = 0. In generic case, the latter possesses precisely ''Jl = {N — 1)^ 
(Bezout's number) zeros in C^. Constructive computation of Vx{g) is possible with the 
aid of an analogue to the division process utilizied in the univariate case. Choose the set 
of power products in X: 

{M,{X)}fJ = {x{'xf\0 < Ji < iV - 1,0 < j2 < 2(iV - ji - 2)} . (12) 
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For instance, one has for N = 5: 



{1, 
xf} 

We will call the reduction of the polynomial Aie{X)g{X) modulo dg/dxi and dg/dx2 its 
representation in the form 

Mi{X)g{X) = beoMo{X) + ... + be,<n-iM<n-i{X) (14) 
+ qii{X)dg/dxi + qe2{X)dg/dx2, 

with {q£i{X),q£2{X)} C C[X]. Theoretical possibility of such a representation as well as 
constructive algorithms for its implementation are discussed in [Ij. We note just only that 
in case of reducibility, the coefficients bij can be expressed as rational functions of the 
coefficients of g{X). Reorder the set f|T2l) in such a manner that Aio = 1, A^i = xi,Ai2 = 
X2 and make the matrix from the coefficients of the reductions f[T^ for i G {0, . . . , 9T — 1}, 
i.e. for all the power products from f lT^ : 

®=NJ~=V (15) 
Denote by ^g^j the cofactor to the corresponding entry of the last row of 03. 
Theorem 3 One has 

Vx{g)=det^. 

The polynomial g{X) possesses a multiple zero A = (Ai, Aa) e (i.e. the zero for which 
g = 0, dg/dxi = 0, dg/dx2 = 0) iffdefiB = 0. Under this condition, the multiple zero is 
unique if ^<yi,i 7^ 0; in this case it can be expressed as 

Ai = ^?n,2/23^,i, A2 = 03^n,3/«ot,i- (16) 



X2, 
X1X2, 
X\X2, 



■^21 
2 

X1X2, 
1-^25 



■^21 



4 
4 

XiX2^ 



X 



2' 



X 



21 



(13) 



Schur formula. Subsequently we will frequently use the following Schur complement 
formula for the determinant of a block matrix [5] : 

det g ^ = det U det (T - SU"^ V) , (17) 

here U and T are square matrices and U is non-singular. 

3 Distance between a quadric and a linear surface 

We treat the equations of the surfaces in the form ([1]) and ([2]) and assume the columns 
Ci, . . . , Cfc to be linearly independent (the latter results in the restriction k < n). Compose 

the matrices C [Ci, . . . , Ck] and 

G =^ C^C, (18) 

i.e. G is the Gram matrix for the columns Ci, . . . ,Ck- Due to imposed restriction on 
Ci, . . . ,Ck, the matrix G is nonsingular. 
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Theorem 4 The condition 



< 



Ai Bi C 

Bf -1 O 

o o 



-1)^ ^, if Ai is positive definite, 
-1)", if Ai is negative definite 



(19) 



is the necessary and sufficient one for the linear surface ^ to intersect the ellipsoid ([iP; 
in this case one has d = 0. If this intersection condition is not satisfied then the value d"^ 
coincides with the minimal positive zero of the equation 



( 



/i 



v 



Ai 



5i 
-1 + iiz 



C 



I 



(20) 



provided that this zero is not a multiple one. 

Proof. I. Finding the intersection condition. Let us first find the critical value 



oil y{X) = X^AX + 2B^X - 1 in the surface C^X 
Lagrange function 



The critical point of the 



X^AX + 2B'^X - 1 - uiC^X - ... - z/fcCjX 
satisfies the system of equations 

2AX + 2B-C [ui, 



Therefore 



X 



-A-^B + -A-^C[ui, 



C'X 



with 

Substitution of 



[z/i,...,i/fc]^ = 2(C^A-iC) 'C^A-^B . 
into ([21]) yields 

Xe = -A-^B + A^C {C^A-^Cy^ C^A-^B 
and the corresponding critical value of V{X) subject to X = O equals 

V{Xe) = -{B^A-^B + 1 - B^A-^C{C'^A-^C)-^C'^A-^B) . 
With the aid of Schur formula (ITTll one can transform the last expression into 



(21) 
(22) 



V{X, 



C^A-^C C^A-^B 
B^A-'C B^A-'B + l 





ABC 




B^ -1 O 


o o 



det(A) det(C^A-iC) 



det(C^A-iC) det(A)det(C^A-iC) ' ^^^^ 

If V{Xe) = then the linear surface ([2]) is tangent to the ellipsoid ([1]) at X = X^. 
Otherwise let us compare the sign of V{Xe) with the sign of V{X) at infinity. These signs 
will be distinct iff the considered surfaces intersect. If A is positive definite then Voo > 0, 



^To simplify the notation we will type matrices A and B without their indices. 
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det(A) > and det(C'^A-iC) > 0. Therefore, V{Xe) < iff the numerator in ([23D is 
negative. This confirmes (!T9|) . The case of negative definite matrix A is treated similarly. 

II. Distance evaluation. Using the Lagrange multipliers method we reduce the 
constrained optimization problem to the following system of algebraic equations 

X -Y - XAX -\B = (24) 

X-Y + ^C[\i,...,XkV = (25) 

X^AX + 2B'^X -1 = (26) 
C^F = O . (27) 

We introduce also a new variable responsible for the critical values of the distance function: 

z-{X - Yf{X -Y)=0 . (28) 

Our aim is to eliminate all the variables from the system (!24|) - (!28|) except for z. We 
express first X and Y from (12^ and (125|) (hereinafter I stands for the identity matrix of 
an appropriate order): 

X = -A-^B-^A-'C[Xi,...AkT (29) 
Y = -A-ifi-^(A-i-AI)C[Ai,...,A,]^. (30) 

Then we substitute ( |30l) into (1271) with the aim to express Ai, . . . , via A. This can be 
performed with the aid of the following matrix 

M Ic^A-^C - C^C = /iC^A-^C - G, (31) 
A 

with G defined by ( ITSll and jji'^= 1/X. Indeed, one has 

M[Ai,...,Afc]^ = -2C^A-i5 (32) 

and, provided that M is non-singular, 

[Ai,...,Afc]^ = -2M-^C^A-^5. (33) 

Now substitute (!33l) into ( l25l) and then the obtained result into (!28l) : 

z - B'^A-^Cyr^GMr^C^A-^B = 0. (34) 

Equation ( |34l) is a rational one with respect to the variables fi and z. 

To find an extra equation for these variables, let us transform (l26l) using (l29l) and ( l33l) 



= X^AX + 2B^X -1 

= -B^A-^B - 1 + fiB^A-^CM'\fiC^A-^C - G + G)M'^C^A-^B. 

Using ( 13T|) and the last equation takes the form 

^(;^, ^) t/ _i + _ B^A-^B + nB^A-^CM-^C^A-^B = 0. (35) 
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Therefore, system fl21]) - fl25]) is reduced to fl51]) - fl5S]) . It can be verified that the left-hand 
side of (IMI) is just the derivative with respect to /i of that of (135!) and, thus, it remains 
to eliminate /i from the system 

^(//,z) = 0, ^'^{fi,z) = 0. 

Taking into account Corollary [T] from Sect. [21 one can perform this with the aid of 
discriminant - and that is the reason for its appearence in the statement of the theorem. 
Schur formula f|T7|) helps once again in representing ^'(/i, z) in the determinantal form: 



^{fi,z) 



A 


B 


C 




-1 + fiZ 


o 




O 






A C 






-G 





A 

B^ 



B 

-1 + ^z 



C 



G 



det(A) det(M) 



(36) 



III. Finding the nearest points on the surfaces. Once the real zero z = z^, oi 
fl20l) is evaluated, one can reverse the elimination scheme from part II of the proof in order 
to find the corresponding points X* and on the surfaces. 

For z = z^, the polynomial in /i standing in the numerator of ( !36ll 



A 

B^ 



B 

-1 + fiZ 



C 



(37) 



has a multiple zero fi = fi^:. Provided that the multiple zero is unique, it can be expressed 
rationally in terms of the coefficients of this polynomial (and consequently in z^,) with 
the aid of (j8]). We substitute this value into (|3T1) then resolve the linear system ( l32l) with 
respect to Ai, . . . , and, finally, substitute the obtained values into (|29l) and ( l30l) . 

However, this algorithm fails if for fi = fi^, the matrix M becomes singular. For 
explanation of the geometrical reason, one may recall that the distance between the 
surfaces may be attained not in a unique pair of points. 

We avoid this case by imposing the simplicity restriction for the minimal zero of J^{z) 
in the statement of the theorem. 

IV. Nonsingularity of the matrix M. In accordance with Theorem |51 the polyno- 
mial J^{z), being the discriminant of $(yU, z), permits the linear representation 



T{z) = v{fi, z)^ + uifi, z)(!>'. 



(38) 



with the polynomials {f (//, z),u{fi, z)} C M[/i, z] satisfying the degree restrictions: deg^ u < 
deg^$,deg^t; < deg^ 

If z = z^: stands for the zero of J^{z), then $(//, z^) and $^(/i, z^,) possesses a common 
zero fi = fj,:^. Differentiate fl38|) with respect to z: 



r(z) 



and substitute n = fi^, z 



r(z. 



(39) 
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We intend to prove that z*) = 0. For this aim, differentiate f l38p with respect to /i: 



and substitute /i = /i^,, z 



= 



52$ 



or 



52$ 



(40) 
(41) 



The second alternative from f HT]) has the meaning that the zero /i = /x* is of multiphcity 
/c greater than 2 for $(/i, z*). In this case, one has from 



(42) 



= i;(/i, 2:*)$(/i, z*) + M(/i, 2:*)$^(/i, z^) ^ 
u{fi,z^)^'^{fi,z^) = -v{fi,z^)^{fi,z^). 



Since the multiphcity of /i = /i* for $^(/i,2;*) equals A; — 1 it follows from fH2|) that 
its left-hand side is divisible by (yU — yU*)'^ while one of its factors is divisible at most 
by (/i — n^Y^^. Consequently, u{ij,,z^) is divisible by /i — /i* and hence u^n^.z^) = 0. 
Therefore, in any case, the condition (HUj) implies that u{fi^,z^) = 0. Formula ( 15^ yields 
then that ^-"'(2;*) = f (/i*, -2*)9<l'/c?2|(.^^ and provided that z^ is a simple zero for J^{z), 
one has J^'{z^) 7^ which results in d^/dz\(^^^ 7^ 0. To obtain the expression for the 
last derivative, let us differentiate the determinantal representation (j3j 



9$ 

dz 



A 



B 



C 



det(A) det ( --G - C' A-'C 



T 



c 

o 
1 



G 



fc+i 



A 



C 



--G 



1)^'+V det(A) det(G + /iC'^'A-^C 



-1 



/idet(A) det(M). 



Since d^/dz for fi = z = z^:, the matrix M should be nonsingular for these values. 
□ 

Corollary 3 // the system of columns Ci, . . . ,Ck is orthonormal then, by transforming 
the determinant in [2D\) . one can diminish its order: the expression under discriminant 
can he reduced into 

Ai - /iCC^ fii 

fif -l+liz 

Corollary 4 Let if G M'^ he the given column. The condition 



(43) 



< 



Ai 




C 






-1 




X 


C^ 


-H 








-1)'^ ^, if Ai is positive definite, 
-1)", if Ai is negative definite 



is the necessary and sufficient one for the ellipsoid (Qj) to intersect the affine suhspace 
C'^X = H . If this condition is not fulfilled then the square of the distance hetween the 
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ellipsoid and the linear manifold equals the minimal positive zero of the polynomial 

I 



\ 



Ai 



Bi C \ 
-H Ig . 



(44) 



provided that this zero is not multiple one. 



□ 



Proof is similar to that of Theorem HJ 

Example. Find the distance to the Xi-axis from the elhpsoid 

7 + 6 X2 + 5 Xg — 4 X\X2 — 4 0:2X3 — 37 xi — 12 X2 + 3 X3 + 54 = . 

Solution. One may choose here Ci = [0,1,0]-^, C2 = [0,0, 1]-'", then the determinant 
takes the form 



-7/54 1/27 

1/27 -1/9 -fi 

1/27 

37/108 1/9 



37/108 

1/27 1/9 

-5/54 - fi -1/36 

-1/36 -l+i^z 



Equation ([20]) 

T{z) 



= 516019098077413632 2^ - 15034745857812486912 z'^ 
+ 95300876926947983328 z^ - 421036780846089455856 z 
+ 237447832908365535785 = 



has two real zeros: Zi ^ 0.05712805 and Z2 ~ 22.54560673. Hence, the distance equals 
/^T ^ 0.23901475. 



Corollary 5 The square of the distance from the point Xq to the ellipsoid (J\) coincides 
with the minimal positive zero of the equation 



T{z) \ det 



Ai 5i 

Bl -1 



+ /i 







(45) 



-I ys.0 
X(f z — X^Xo 

provided that this zero is not a multiple one and X^AiXq + 2BJXq — 1 7^ 0. 

The square of the distance from the origin X = O to the ellipsoid (Q]) coincides with 
the minimal positive zero of the equation 

J^(z) =^ {{fiz - 1) det(Ai - fil) - 5i^adj(Ai - /il)5i) = (46) 
provided that this zero is not a multiple one. Here adj stands for the adjoint matrix. 

Remark. For large n, one can compute det(Ai — /il) and adj(Ai — yul) simultaneously 
with the aid of the Leverrier-Faddeev method [3]. 

Remark. For the case Bi = O, one gets J^{z) = V{f) [z"'f{l/ z)]^ with /(/i) = 
det(Ai — /il). This corresponds to the well-known result that the distance to the el- 
lipsoid X^AiX = 1 from its center coincides with the square root of the reciprocal of the 
largest eigenvalue of the matrix Ai. 
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We exploit the result of the last corollary to elucidate the importance of the simplicity 
restriction imposed on the minimal positive zero for T{z)\ this assumption will also appear 
in the foregoing results. 

Example. Find the polynomial (H5|) for the ellipse /A + y"^ = 1 and the point (xq, yo)- 
Solution. The polynomial J-'{z) from for the ellipse x^/a^ + y^/fe^ = 1 is computed 

as 



J-(z) = V^(^fi^-{a^ + b^-xl-yl + z}fi' 



+ i-a%'(^ + ^-l]+zia' + b')\fi-a%'z 



62 



which for our particular case a = 2, 6 = 1 yields (up to a factor 1/256) 

T{z) = 9z^ -Q{2xl + 7y'^ + 15) z^ 

+ (-2 x^Q + 73y^ + 62 x^^ - 90 a;^ + 270 yl + 297) 
+ (-56 - 360 - 62 - 248 + 4 + 270 xl 
- 90 xly^ - 30 x^^yl + 140 xly^ - 360) z 

+ A{xt + 2xlyl + y^,-6xl + Qyl + 9){xl/A + yl-l)' . (47) 

Let us evaluate its zeros for yo = 0, i.e. for the points in x— axis: 

J^{z) = {z- (xo - 2)2)(z - (xo + 2)2)(3 z-{3- x^))l 

Multiple zero Z2 = I — Xq/3 is positive for xq G [0, a/3 [. Moreover, for these values of Xq, 
zero Z2 is the minimal one for J-'{z). Nevertheless, for xq > 3/2, the square of the distance 
from (xq, 0) to the ellipse equals zi = (xq — 2)^. Explanation of this phenomenon is as 
follows: the multiple zero Z2 corresponds to the pair of points (4xo/3, ±^^1 — 4xg/9) in 
the ellipse. These points are real for xq G [0, 3/2 [ and imaginary (complex-conjugate) for 
Xo > 3/2. 

To conclude this example, let us illuminate the relationship of the stated metric prob- 
lem to an ancient one concerning conic sections. Let us estimate the number of real 
zeros of the polynomial (H71) . For this purpose, the sign of discriminant of polynomial is 
significant. One has 

Drawn in the (xq, |/o)— plane, the curve \E'(xo,2/o) = consists of three branches: the 
coordinate axes and the curve known as astroid (marked in red in Fig. [1]). The latter 
was first treated by ApoUonius in the 3rd century BC, in connection with the problem 
of finding the number of normals drawn from the given point to the ellipse. In terms 
of the zeros of polynomial P7|) . the solution is as follows: for the points (xo,?/o) inside 
the astroid the polynomial J-'{z) posseses four real zeros, for those outside - two. The 
exceptional points lie in the axes: one gets four real zeros for corresponding J^{z) (with 
two of them becoming negative outside astroid). 

To complete the present section, we provide estimations for degrees of polynomials 
J-'{z) appeared in the above results. 



Theorem 5 For the polynomial from /[2^) . one has generically deg J-" = 2k. 
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Figure 1: 



Proof. For simplicity, we will treat the case where the columns Ci,...,Ck are or- 
thonormal. We expand first the polynomial under the discriminant sign in powers of 

_ _ ^ Ai 5i C 



k+l 



Ai C 



-1 



l//il 





Ai 


C 






l//il 



Here I stands for the identity matrix of order k. The leading term of J^{z) = z)) 
coincides with 

In order to evaluate the degree of the last expression w.r.t. variable z, we may exploit the 
formula ([9]). For this aim, it is necessary to find the degree of the polynomial under the 
discriminant sign w.r.t. fi. Application of Schur formula (ITTj) in the way corresponding 
to (03]) yields 



ZfX 



k+l 



Ai C 
l//il 



2;/idet(Ai - yuCC^ 



which is not useful for our purpose since the matrix CC"^ is singular ii k < n. Let us use 
Schur formula in an alternative way: 



= zfidet Ai det(I-/iC' A^^C). 

The last determinant is of the order k with all of its entries depending linearly on /x. We 
expand it in decreasing powers of /i: 

= det Ai det(C'^A7^C) + . . . 

Since, by the assumption, matrix Ai from the equation ([1]) provides an ellipsoid, all the 
matrices Ai, A^^ and C^^A^^C are sign-definite. Therefore, their determinants do not 
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vanish and the leading term of T{z) equals generically 

z^^{det A^f^V^{deti\ - /iC^A^^C)) = z^^{dei AifV^{det{Ai - /iCC^)). 

□ 

Corollary 6 For the polynomial from the leading term equals generically to 

^^"(det Ai)2l)^(det(Ai - /il)). 

4 Distance between quadrics 

Consider first the case of surfaces ([1]) and ([3]) centered at the origin: Bi = O, 82 = O. 

Theorem 6 The surfaces X'^ AiX = 1 and X"^ A2X = 1 intersect iff the matrix Ai — A2 
is not sign- definite. If this condition is not satisfied then the value d"^ coincides with the 
minimal positive zero of the equation 

J^{z) P;,(det(AAi + {z- A)A2 - \{z - A)AiA2)) = (48) 
provided that this zero is not a multiple one. 

Proof. I. The intersection condition can be found as an exercise in the problem book 
[7]. We just repeat here the arguments. 

Since the equation X^AiX = 1 provides an ellipsoid, the matrix Ai is positively 
definite. 

Let there exist a point X = Xq G such that X^f AiXq = 1 and X^f A2X0 = 1. Thus 
Xq{Ai - A2)Xo = for Xq ^ O. Therefore, Ai - A2 is not sign-definite. 

Conversely, if Ai — A2 is not sign-definite then there exists Xq 7^ O such that X^f (Ai — 
A2)Xo = or, alternatively, X^AiXq = X^f A2X0. Multiply the latter by a scalar with 
t e R: t^X^AiXo = t^Xff A2X0. Set t = I/a/xJA^Xq (the radicand is positive due 
to the positive definiteness of Ai). The point X = tXg is an intersection point of both 
manifolds since 

X^AiX = t^xjAiXo = 1 and X^A2X = t'^X^A2Xo = t^X^AiX^ = 1. 

II. If the intersection condition is not valid, then the distance problem becomes non- 
trivial and we apply the Lagrange multipliers method for the objective function in the 
form 

(X - Yf{X -Y)- Ai(X^AiX - 1) - A2(r'"A2r - 1). 
The corresponding system of algebraic equations is as follows 

X -Y - AiAiX = O, -X + F - A2A2F = O, (49) 

X^AiX = 1, Y^A2Y = 1. (50) 

This system yields 

(A1A2A2A1 - AiAi - A2A2)X = O, (51) 

(A1A2A1A2 - AiAi - A2A2)F = O, (52) 
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and 

AiAiX + A2A2F = O, (53) 
X-Y = AiAiX. (54) 

Matrices of the systems flSTl) and fl52l) differs only by transposition, and therefore the 
determinants of these matrices are equaL Their common value should be just due to 
the fact that we are looking for nontrivial solutions of homogeneous systems: 

det(AiA2AiA2 - AiAi - A2A2) = 0. (55) 

Let us introduce the matrix 

M I _ Ia^^ - ^A,' (56) 
Ai A2 ' 

and the vector 

Z t/ X - Y. (57) 
Using this notation, the equations (jSTj) and ( 152|) can be rewritten into equivalent form 



(58) 



while the conditions fl50|) in the form 

1 



^Z'^P^jZ = 1 for J G {1,2}. (59) 



Let us introduce a new variable z responsible for the critical values of the distance function 

z = {X~YY{X-Y)^ Z^Z (60) 

^Z^A-,'Z + ^Z^A,'Z # Ai + A2. 
Ai ' A2 ' 

Thus, we have eliminated the variables X and Y from the system fl49p -f l5Up with the 
resulting equations assuming the form ( 155|) and ( 160|) . To deduce an extra equation, one 
should start with the identity 

M-adj(M) =I-detM. 

By differentiation this as to Aj, one obtains 

^^5adj(M) 9detM^ 

Multiply this by from the left-hand side and by Z from the right-hand side, with Z 
standing for any nontrivial solution to the system f l58|) : 

T^dm , r 5adj(M) (9detM ^ , , 

Taking into account (I58p and symmetry of the matrix M, one arrives at 

z^m = {MZf = O, 
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and therefore identity fl6ip turns to 
or, in view of 

^Z^Aj'^dm)Z = ^^Z^Z. (63) 
Now, our aim is to prove that 

adj(M)Z = 7Z (64) 

for a certain scalar 7. Indeed, 

adj(M)MZ = O ^ M(adj(M)Z) = O. 

If rank(M) = n—1 then any solution U to the system of homogeneous equations Mf/ = O 
should be equal just a multiple of Z; therefore 

adj(M)Z = 7^. 

The case rank(M) < n — 1 is trivial since adj(M) = 0„xn- (It can be proved that in any 
case 7 = Mil + M22 + . . . + M„„ with Mjj standing for the cofactor to the corresponding 
entry of M.) 

Hence, the formula ( 163|) is transformed into 



wherefrom one can deduce (with the aid of (l59i) ) that 

9detM _ adetM 



(65) 



Recalling now that Ai and A2 are connected via condition (!60|) . we substitute Ai = 
2; — A2 into ( 165|) and obtain 

adetM dA2 _ SdetM ^ ddetM_^ 
8X2 dXi 8X2 8X2 

Thus, the process of elimination of variables from the system fl50p - (|S2]) and (jSO]) terminates 
when we get the two equations: the first one is 

det (l - - ^A,') = 0, (66) 



z-Xo ' X 



2 



while the second is obtained from this by differentiation its left-hand side as to A2. Elim- 
ination of A2 from these equations can be performed in the traditional manner, i.e. via 
discriminant. Utilizing the result of Corollary [1] from Sect. we turn from the ratio- 
nal functions to polynomial ones. Multiplication of (!66|) by det(AiA2) and substitution 
A = I/A2 completes the proof. 

III. To find the nearest points on the quadrics we suggest the following approach. 
Once the real zero z = oi is evaluated, one can find the corresponding value 
A = A^, which is a multiple zero for 

G(A, z^) =^ det(AAi + (z, - A)A2 - X{z, - A)A2Ai). 
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Under the assumption of the theorem, this zero is unique and can be expressed rationally 
in terms of the coefficients of the polynomial G(A,z*) with the aid of (fT6|) . Futhermore, 
the coordinate column X* of the point on the quadric X-^AiX = 1 is a solution for the 
system of homogeneous equations 



(A*Ai + (z^ - K)A2 - K{z^ - A*)A2Ai)X 



(67) 



which possesses an infinite number of solutions since its determinant vanishes. From the 
solution set one should choose a representative satisfying the condition X^AiX = 1. Due 
to symmetry of the problem, there exists a pair of such solutions. 

Similarly, the coordinate column for the point in the second quadric satisfies the system 

(A,Ai + {z, - K)A2 - K{z, - A,)AiA2)F = O. (68) 

Recall that the matrices of the system f E7|l and f l68p differ only by transposition and 
in order to solve both systems (|67|) and ( 168|) it suffice to treat the rows and the columns 
of the matrix adjoint to the matrix 



A*Ai + (z* — A*)A2 — X*{z^ — A*)A2Ai. 



Indeed, Xj equals just a multiple of any nonzero row of the matrix adj(M*) while 
coincides with a multiple of any nonzero column of adj(M^,). By a suitable selection of the 
mentioned multipliers, one can provide the fulfilment of the conditions X-^AiX = 1 and 
Y^AiY = 1. The obtained pairs of points should be adjusted according to the condition 



(x« — Y^)'^(x^ — y^,) 



Example. Find the distance between the ellipses 



□ 



10 Xi — 12 X1X2 + 8 X2 = 1 and Xi + X1X2 + = 1 • 



Solution. Here 



Ai 



10 -6 

-6 8 



A, 



I 1 



and the matrix Ai — A2 is positive definite. Thus the ellipses do not intersect. 
Compose the determinant from Theorem [61 



G{X,z) 



det(AAi + {z- A)A2 - X{z - A)AiA2 



7X'^ + {-7z + 9)X + z 



X^+{z-^]X + ^z 



149 A 



2z - ^ ] X + ^z 



-2A2 



5X^ + {-5z + 7)X + z 



(69) 



33 A* + -66z + —]X''+\33z 



2 J 



61. + ^ 



27 



45 



-^z^ + \X + -z 



+ 

The discriminant of this polynomial w.r.t. A equals 

3 



Hz) 



16 



2^(936086976 z"" - 10969697376 + 50706209664 z^ 
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-115515184664 + 130176444432 z"^ - 59826725574 z + 2866271785) 
Its positive zeros are as follows: 

2*^ 0.053945666, 1.3340583883, 1.95921364, 2.8785867381. 



Hence, d = ^ 0.23226206. 

To find the nearest points on the given ellipses, establish first the multiple zero of the 
polynomial G(A, z*) with the aid of Theorem [1] 



-725274^5 + 1455894 2^ + 



11286981 , 26486523 , 42000075 



z^ + 



A 



17591706^4 - 109992894^3 + 



450450691 ^ 315606253 77466805' 
z'^ z H 



Substitution z = z^ yields A = A,, ~ —0.13576051. 

Secondly, for the obtained pair of values z* and A* the determinant f l69|) vanishes and 
therefore both systems f l67|) and fl68|) : 

M,x = o, yfy = o 



posses nontrivial solutions. To find these solutions, take the first row and the first column 
of the matrix adj(M^,) 



X 



Y 



2X' 



(2^-f)A 



7X^ + {-7z + 9)X + z 
\^-(z-f)\-lz 
7\'^ + {-7z + 9)\ + z 



z=Zt ,A=A* 



-0.8579069 

-0.9876166 
-0.8836615 
-0.9876166 



Each such point defines a line passing through the origin. To find intersection points with 
the corresponding ellipses, one should make normalization 



X^ 



±x 



VX'^'AiX 



± 



-0.3838312 
-0.4418639 



±Y 



Vy^a^y 



± 



-0.5449964 
-0.6091105 



These formulas provide two pairs of nearest points in the ellipses. 
Let us now treat the general case of manifolds position. 

Theorem 7 The surfaces X^AiX + 25f X - 1 = and X^ A2X + 2SjX - 1 
intersect iff among the real zeros of the equation 



$(2) t/ Vy^ \ det 



B2 

■1-z 



A 



Ai 5i 

BI -1 







there are the values of different signs or 0. If this condition is not fulfilled then the value 
d"^ coincides with the minimal positive zero of the equation 



Hz) 



def 



V 



det ( /i 



Ai B, ' 
_BT -1, 
A2A1 A2B1 
BjAi B^Bi - niHiZ 



+ At2 



A2 B2 
BJ -1 



(70) 



provided that this zero is not a multiple one. 



17 



Proof. We sketch it as it is similar to that of Theorem |H Intersection condition is 
a result of the following considerations. Extrema of the function X'^ K2X + 2B\X — 1 
on the ellipsoid ([1]) are all of the similar sign iff the surfaces ([1]) and (|3]) do not intersect. 
We state the problem of finding the extremal values of ViX^ = X'^A2X + 2i?JX — 1 
subject to ([T]), then apply the Lagrange multipliers method and finally eliminate all the 
variables except for z from the obtained algebraic system coupled with the equation 
X^A2X + 2B^X - 1 - z = 0. 

To prove the second part of the theorem, we take the matrix M defined by fl56|) . while 

and transform the equations of the system (jl]) into 

-A,'B2 - ^A,'M-'Q, (71) 
l = 0forjG{l,2}, (72) 

(73) 

On multiplying equations fl72|) by and using fl73l) . we deduce that 

- XiB^A^'Bi - X2B^A^'B2 - Q^M-'Q - Ai - A2 + 2 = 0. (74) 

It can be verified that the derivative of the left-hand side of fTMl) with respect to Aj 
coincides with that one of fl72l) . Substitution /xi = I/A2, /i2 = 1/Ai and the use of Schur 
formula ( |T7|) enable one to reduce (1741) to the determinantal representation from ( 1701) . □ 
Example. Find the distance between the ellipsoids 

7 + 6 X2 + 5 X3 — 4 X1X2 — 4 X2X3 — 37 xi — 12 X2 + 3 X3 + 54 = 

and 189 x\ + x\ + 189 0:3 + 2 X1X3 — X2X3 — 27 = 

and establish the coordinates of their nearest points. 

Solution. Intersection condition from Theorem [7] is not satisfied: the 6th degree poly- 
nomial ^{z) has all its real zeros positive. To compute the discriminant we use the result 
of Theorem [3] with the matrix of the order 16 constructed for the set f|T3|) . Its deter- 
minant is the 24th degree polynomial J^{z) with integer coefficients of the orders up to 
10^*^®. It has eight positive zeros 

zi ^ 1.3537785, Z2 ^ 3.5509348, . . . , ~ 111.7480312. 

Thus, the distance between the given ellipsoids equals ^/zl ^ 1.1635198. 

For the obtained value of zi, polynomial in /ii and /i2 from ( [70|) possesses a multiple 
zero which can be expressed rationally in terms of zi with the aid of the minors of 03 via 
formulas ( fT6l) . Substitution of the obtained values Ai ~ 5.75593612, A2 ~ —0.45858332 
into (1711) yields the coordinates of the nearest points on the given ellipsoids: 

X ^ [1.5203947, 1.5098600, 0.1262343]^, 

Y ^ [0.3610045, 1.4849072, 0.0315226]^. 



X = -A^'Bi + — A^^M-iQ, Y = 
Ai 

-BjAj'B, + Ig^M-^ATiM-ig - 

j 

z - Q^M^^Q = 0. 
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5 Parameter dependent surfaces 



The problem of distance estimation between moving objects in 3D space is of importance 
to astronomy, robotics and computer graphics. To illuminate the perspectives of the 
approach developed in the previous sections for such problems dealing with quadrics, we 
will treat the following problem. 

Find the distance from the point Xq G M" to the nearest point of the family of ellipsoids 
in M" 

^^X^Ai{t)X + 2Bf{t)X -1 = tG[a,6]} (75) 
with the coefficients of Ai(t) and Bi{t) polynomially dependent on the parameter t. 



Theorem 8 The square of the distance from Xq to ( [73] ] coincides with the minimal pos- 
itive zero of one of the equations 

^{z) Vt{J^{z, t)) = 0, T{z, a) = 0, J^{z, h) = 0. 
Here J^{z^ t) is a polynomial and the mentioned zero is not a multiple one. 

In short: the stated problem can be solved with the aid of iterated discriminant. 
Proof. For any given value of t, the square of the distance from Xq to the correspond- 
ing ellipsoid of the family fl75|) is evaluated as a zero of the equation fH5l) 



J^(2,t) = 0. (76) 

Due to imposed restrictions on the coefficients of the family, J-" is a polynomial function 
in t. Equation fl76l) can be treated as defining an implicit function z{t). It is known 
that zeros of a polynomial are continuously different iable functions of the coefficients of 
this polynomial (except for the coefficient specializations annihilating the discriminant) 
0. Consequently, for any zero z = z^{t) of (1751) there exists the derivative dz^(t)/dt. 
Differentiation of the equality J-'{z^,{t),t) = with respect to t results in 

dz,{t) d^^r. 
dz' dt ^ dt " ' ^ ' 

here the partial derivatives are evaluated at z = z^:{t). 

For t G [a,b], the minimum of the function z*(t) is attained either at the end points 
of the interval or in the stationary point t = t at which dz^^/dt = 0. In the latter case, it 
follows from fl77|l that 

(-) 

at t = t. The two conditions fl76|) and fl78l) provide an algebraic system with respect to 
both variables z and t. One can eliminate the variable t with the aid of discriminant. □ 
Example. Find the distance from the point (—10, 10) to the family of ellipses 



{x-ty , (y-t{t-A)Y 



t G 



+ 16 

Solution. We skip the expression for J'{z,t). 

Vt{J^{z,t)) = z^(3z+ 16888)2(9 2^ -4080^ + 333376)2 

X ( 16777216 _ 24039653376 + 15135396003840 
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- 5551772745220096 + 1322366761276505856 / 

- 215049198876048266976 + 24423380307243182292153 

- 1952292050779441220868024 z^ + 109783307459960901970173936 / 

- 4304075084512715479517135104^3 

+ 113714594973157300688449668864 z'^ 

- 1830069428535779484150176987136 z 
+ 14265422520155306699255826485248)3 

X (2304 / - 3774720 z'^ + 2645308000 / - 1058624029488 z^ 

+ 266900597798217 z^ - 42785419475837458 z^ 

+ 4100511694812810849 - 202905147887926860744 z 

+ 3648597980765724103824). 

The minimal positive zero of the last factor is 2* ~ 37.70933565. The distance to the 
family equals ^/z^ ^ 6.140792755. One can find the ellipse of the family at which the 
distance is attained via the traditional application of the multiple zero evaluation formula 
(ED: U ^ -1.9680233599. 

6 Conclusions 

We have treated the problem of distance evaluation between algebraic surfaces in via 
inversion of the traditional approach: 

nearest points — )■ distance . 

This has been performed via introduction of an extra variable responsible for the critical 
values of distance function and application of Elimination Theory methods. Such an 
approach was first suggested in [10] for the general polynomial optimization problem 
in M". Its employment for the distance evaluation problem for quadrics has led to the 
result which happened to be surprisingly unpredictable for us: the discriminant is fully 
responsible for everything. With its help it is not only possible to deduce a univariate 
polynomial equation for the square of the distance but also to express (Theorem [7]) the 
necessary and sufficient condition for the intersection of the surfaces. 

The major advantage of this approach over the traditional scheme is that the problem 
is reduced to evaluation of a single zero of a univariate algebraic equation instead of 
dealing with multidimensional constrained optimization problem. Moreover, introduction 
of an extra (distance) variable z into the problem provides one with a nice (i.e. rational) 
parameterization of the nearest points coordinates. 

Several problems have remained for further investigation, among them estimation of 
the degree of polynomial J^{z) constructed for the problems of Sect. HI The conjecture 
is that deg J^{z) = n{n + 1) for J^{z) from fHHj) and that deg ^-"(2;) = 2n{n + 1) for F{z) 
from ( iTOj) : with these estimations valid on excluding some extraneous factor (e.g. in case 
of dlH]), this factor is just z'"("-i)). 

The proposed approach might be especially useful for the optimization problems con- 
nected with the parameter dependent surfaces like the one treated in Sect. [5] or for the 
multidimensional pattern recognition analysis. 
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